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Abstract 

We study the energy conditions and geodesic deformations in Bertrand 
space-times. We show that these can be thought of as interesting physical 
space-times in certain regions of the underlying parameter space, where 
the weak and strong energy conditions hold. We further compute the 
ESR parameters and analyze them numerically. The focusing of radial 
time-like and radial null geodesies is shown explicitly, which verifies the 
Raychaudhuri equation. 
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1 Introduction 



The Schwarzschild metric, discovered nearly a century ago, remains one of the 
simplest yet most profound solutions of Einstein's field equations. One of the rea- 
sons for the popularity of the Schwarzschild solution among relativists is that it 
provides a realistic scenario to describe closed, stable orbits of planets and other 
heavenly objects. It is however well known that there are other solutions of Ein- 
stein's equations which can also describe such stable, periodic motion. One class 
of examples was discovered in a remarkable paper by Perlick and these were 
named "Bertrand space-times" (BSTs). Indeed, Perlick's classification general- 
izes the well known Bertrand's theorem [2] in Newtonian mechanics (an excellent 
exposition can be found in [3]) to general relativity. The former theorem states 
that the harmonic oscillator and Kepler potentials are the only spherically sym- 
metric potentials for which bounded orbits are periodic. Based on the standard 
deductions of the Bertrand's theorem there was an attempt to generalize its form 
in the special relativistic case [4J. Perlick's work determines all static, spherically 
symmetric space-times in the most general case where one can have stable, closed 
orbits. 

Apart from being interesting from a purely theoretical perspective, BSTs 
might also be relevant for other reasons. For example, one possibility may be 
to model a realistic space-time that allows for bounded, periodic orbits. Al- 
though such a possibility was ruled out in Perlick's original work due to the fact 
that asymptotically flat BSTs (relevant for the motion of objects around an iso- 
lated mass) do not seem to satisfy the weak energy condition (WEC) at infinity, 
asymptotically non-flat BSTs are equally interesting objects, as alternatives to 
black hole space-times. 

This paper studies a class of BSTs from the perspective of the energy con- 
ditions and geodesic deformations. We analyze these aspects and find that in 
certain regions of the parameter space, BSTs do obey the strong and weak en- 
ergy conditions. We further study the Raychaudhuri equation for BSTs and 
confirm the geodesic focusing theorem. 

The paper is organized as follows. In the next section, we briefly review 
BSTs and analyze the energy conditions therein. In section 3, we study radial 
and circular geodesic flows in a class of BSTs and analyze the focusing theorem. 
Section 4 ends with our conclusions and directions for further study. 
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2 Energy Conditions in Bertrand space-times 



Formally, the definition of a Bertrand space-time P, [5] arises via a static, spher- 
ically symmetric Lorentzian manifold (M, g) whose domain is diffeomorphic to a 
product manifold (ri , x S*^ x M with the metric g given by 

where r ranges in the open interval (ri , r2), 5 and are co-ordinates on the 
two-sphere. A and u are some unspecified functions of r to start with. Such 
a Lorentzian manifold is called a BST provided there is a circular trajectory 
passing through each point in the interval (ri , r2), which is stable under small 
perturbations of the initial conditions. 

Starting from this definition, Perlick ^ deduced that there can be two cate- 
gories of BSTs given by: 

ds^ = 



2[1 - Dr^ ± J(l - Dr^y - Kr^ , , 9 . , 9x 

+ m-Dr-r--Kr*\ ^ *^ + ''V"^ + » . (2) 



which will be called the Type I and Type II forms of the BST respectively. The 
parameters D, G and K are real, and /3 must be a positive rational number. 

For mathematical simplicity, we will consider BSTs of type II, defined by the 
metric of eq.(l3]). We wish to first understand what type of matter distribution 
can cause this metric, assuming the Einstein equations to hold. To this end, 
we construct the Ricci scalar and the energy momentum tensor. The general 
expressions are too lengthy to reproduce here, and we will frame our arguments 
based on special cases. Let us begin with the case = for which the type II 
metric of eq.([3]) reduces to 

ds' = --^ + ^ + r\de' + sin' e d<p') , (4) 

where we will take G > to ensure a Lorentzian metric. The Ricci scalar can be 
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calculated to be 







^ ^ 4(1 + Grf - (3^ {AGr (2 + Gr) + 7) 

2r2 (1 + Gr)^ ' 



this diverges at r — )■ and vanishes as r — )■ oo. The stress-energy tensor is 
proportional to 

^ (Gr + 1) /32((Gr + 2)/32-Gr-l) /^^(l - 2Gr) - 2Gr) \ 

^""^V ' r2(Gr+l) ' Ar\Gr + lY' Ar\Gr + lY) 

(6) 

To analyze the energy conditions, it is convenient to introduce an orthonormal 
frame that satisfies 

£//.^e^e^ = (7) 

where 7]^^ = diag(— 1, 1, 1, 1) is the fiat Lorentzian metric. Since the metric of 
eq.(j4]) is diagonal, it is easy to see that a choice of the orthonormal basis is given 
by eit = diag ( -r^ — , ] whence the energy momentum tensor can 

" ^ V ^900 ' v^n V^22 V933 J 

be written as 

T^'' = pe^e^ + pie^er + p^e^e^ + pse^e^ (8) 
and the energy density p and the principal pressures Pi, = 1 ■ ■ ■ 3 are 



1-/3^ _ I3^{2 + Gr) - (1 + Gr) _ _ (3^ {I - 2Gr) 
-^^P'- r2 (1 + Gr) ' P'-P^- 4^2 (i + Gr)' 



It is seen that /3 > 1 is ruled out on physical grounds. /3 = 1 is somewhat 
unphysical, as it implies a vanishing energy density in the presence of non zero 
pressures. For /3 — )■ 1~, the weak energy condition [S], [7], p > 0, p + j9j > 0, 
i = 1, ■ ■ ■ ,3 is satisfied for r < j^. The WEC provides an interesting upper 
bound on r, and G has to be a small positive number for a physically meaningful 
solution in this case. For /3 < 1, the WEC is satisfied for certain intervals of r, 
depending on the choice of G. Specifically, it can be checked that for positive 
values of G (necessary to retain the Lorentzian nature of the metric of eq.(|l]) at 
large values of r), the WEC is satisfied for all r. 

Let us now turn our attention to non-zero values of K, where the situation is 
more complicated. With K 0, the Ricci scalar diverges at r — )■ 0, and in the 
limit r — )■ oo, Roo = —QK(5'^. The energy density and the principal pressures can 



^Here and in the rest of this section, we set 6* = ^ in the final expressions, without loss of 
generality. 
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P = -2 (10) 



be found by introducing an orthonormal frame analogous to the case = 0, and 
we find that 

1 - (3irr2 + 1) 

For r ^ 1, this imphes that the energy density is negative for positive values of 
K. The situation might be remedied by choosing a negative value of K, but note 
that this necessitates, from eq.([3]) that for K = —k where k is a positive real 
number, we require r < 1/ a/k. We can thus choose k <^ 1 so that the positivity 
of the energy density of space-time of eq.(l3]) is guaranteed for a large range of r. 
The analysis of the WEC is similar to the case K = considered earlier. We will 
omit the algebraic details here and simply state the result that setting (3 = 1 for 
simplicity, for a given choice of k (in accordance with the discussion above), the 
WEC is always satisfied for r < 1/ y/R,. 

Before we end this section, we will briefly comment on the strong energy 
condition [6], [7] that follows from eq.([8]) : p + YliVi ^0, p + Pi > 0. We find 
that for the metric of eq.(l3]), 

p+yp^ = (11) 

so that the SEC is satisfied whenever r < 1/ for positive values of G. This 
will be important for us in the next section. 

To summarize, in this section we have studied the energy conditions of Bertrand 
space-times of type II, given by the metric of eq.(l3]). An entirely similar anal- 
ysis can be carried out for the Type I metric of eq.(|2]), although the algebraic 
expressions are complicated. We now proceed to study geodesic flows in BSTs. 



3 Geodesic Flows in Bertrand Space-times 

The kinematics of geodesic congruence in any space-time can be specified by 
three quantities: the isotropic expansion, the shear, and the rotation variables. 
In totality these are generally called the ESR variables, and the evolution of these 
are guided by the Raychaudhuri equations [8]. Treating the geodesic congruence 
as a deformable fluid, one can write the evolution equation of the vector between 
two fluid points. This vector may get deformed as the geodesies flow, and conse- 
quently the vector is called the deformation vector. The Raychaudhuri equations 
connect the evolution of the deformation vector with the curvature of space-time. 
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In general the deformation vector is called ^'^ and its rate of change with respect 
to an affine parameter is given as 



e = 5'^.r, (12) 

where the affine parameter interval in which the rate is measured is supposed to 
be small. Here is a second rank tensor characterizing the time evolution of 
the deformation vector, 

5^ = V.^i^ , (13) 

where is a tangent vector field which serves as the first integral of the geodesic 
equations. Specifically, u"VuU^ = 0, and choosing a suitable affine parameter one 
can make u^u^ = — 1 for time-like geodesies while u^u^ = for a null geodesic. 
Differentiating the expression in Eq. (fT2|) with respect to the affine parameter 
one gets 

= {13^, + B^B\)C ■ (14) 

The Raychaudhuri equations are obtained by writing = —R^^^^^u'^u^S,'^ and 
equating this with eq ffT^ . 

In n space-time dimensions, the general form of the second rank tensor 
can be decomposed into irreducible parts as [7] 

Bf,y = —^Qh^^ + cr^^ + uj^y , (15) 
n — 1 

where h^y = g^y -\- u^Uy for time-like, and 9 is the expansion variable, a^y 
is associated with shear and w^j^ signifies rotation. The physical significance of 
these variables are nicely explained in [7]. One can explicitly write 

= (16) 

a^,y = -{B^y + By^) -Qhf,y, (17) 

2 n — 1 

^t^y = -{B^y — By^). (18) 

and the ESR variables are generally denoted by 0, cr^ and u"^. From Eq. (1131) it 
is seen that if one knows the form of one can calculate B^y, and hence the 
ESR parameters. These are expected to give us information about geodesic fiows 
and their properties in BSTs. 
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3.1 Geodesies in BSTs of Type II : General Considera- 
tions 

We now focus on Type II BSTs, and further simplify the situation by choosing 
6' = 71/2, so that we are on the equatorial plane. In that case we have the metric 

The geodesic equations can now be written down. The first one is obvious from 
the form of the above line element, 

C, (20) 



G + VK + r-2 

where C is a constant of integration. This equation can also be written as 

i+ =0. (21) 

r'^y/l + Kr\G + VK + r-^) 

The other geodesic equations are: 
and 



+ -r(/) = O. (23) 
r 

On a radial geodesic Eq. (120|) holds and more over = 0. For a time-like geodesic 
{u^u^ = — 1), one can write 

4-2 -2 

+ rV = -l- (24) 



G + ^K + r-2 /32(l + irr2) 

Using the above equation one can calculate the value of u*" = r = c/r/dA on a 
radial geodesic. Here A is an affine parameter. The value of u'' comes out as 



yj- = ^ = + Kr2) [^2 [^jK + r-^ + G^ - (25) 

The above equation is for the outgoing radial geodesic directed away from the 
origin. Note that, assuming Kr^ + 1 > (see discussion in section 2), this 
implies that there is a turning point of the outgoing radial time-like geodesies. 
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for [VK + + G) = 1. This implies that there is a maximum value of r 
at which outgoing radial geodesies stop. This is analogous to the case of the 
Schwarzschild black hole, where it is known that such turning points occur for 
non-marginally bound radial time-like geodesies. The other components of the 
tangent vector on the radial geodesic are 

„..|.c(g.VFTF^), .^^1^0. (26) 

Having calculated the relevant components of the tangent vectors u'^ on a radial 
time-like geodesic of Type II BSTs, one can compute the components of the B^^^ 
tensor for radial time-like geodesies. 

For future reference, let us also list the components of the tangent vector for 
the radial null and the circular time-like geodesies. For the former, we obtain 



c(g + VK + r-2 



= (3C^{Kr^ + 1)(g + VK + 
u'^ = 0, (27) 

where C is defined in Eq. f l2U]) . It is interesting to note that for null radial 
geodesic case there is no turning point (for G > as is always assumed in this 
article) as was present for the time-like radial geodesies. This implies that light 
propagating away along the radial direction in BST of type II is not bound to 
return after travelling a finite distance. In the later part of this article we will 
see that although outgoing null radial geodesies do not have a turning point, an 
outgoing radial null geodesic congruence does focus away from the origin. 
For the circular time-like geodesies, we find 



V 

u 



V2r 



VK + r-2 {G + VK + r-2) 



u''' = W — ? — 7 , ^ ^ (28) 

^ ' " [2r {GVKr^ + 1 + Kr) + l] ^ ^ 

From the above expressions of the tangent vectors one can see that there exists an 
upper bound on the radial distance up to which BST of type II can accommodate 
time-like circular geodesies. The upper limit is given by the inequality 

2r (cVKr"^ + 1 + Kr^ + 1 > . (29) 



7 



If = the above inequality is satisfied for all r. On the other hand if G = 
and K = —k where k > 0, the upper bound is given by 1/a/2k. In general when 
G > and K < Q the upper bound on r has to be evaluated by solving the 
inequality in eq. (l29l) . If both G and K are greater than zero the inequality in 
eq. fl29|) is trivially satisfied but as we have seen in section [2] that in this case the 
WEC and SEC are violated for r > 1. 

To summarize, in this subsection, we have considered the BST of type II 
(eq.©), and calculated the first integrals of the geodesic equation for radial 
time-like, radial null and circular time-like vectors. These can be used in a 
straightforward manner to evaluate the ESR parameters for BSTs, which we 
now turn to. 



3.2 The ESR variables for Type II BSTs 

In this subsection, we compute the ESR parameter for BSTs of type II. Consider 
first the radial time-like geodesies. We start from eq. (fT3ll . from which we can 
write 

dx" 

The non-zero components of i?^^, required to evaluate the ESR variables for the 
radial time-like geodesic flow, in the equatorial plane for the Type II Bertrand 
space-time, are listed below: 



B''. = i^ + V, (30) 



B\ = . , B' 



2r2(G + VK + r-2) 2rVl + Kr^ 



B\ = — _ . B\. 



(3 



B\ = ^^(1 + Kr^)[C^{G + + r-^) - 1] . (31) 
Now from eq.f ll6l) . we get 

/3Vl + Kr'^iGr + ^/TTK^) 



e 



(32) 



2r\ G^{G + + r-2) - 1 



The shear coefficient squared, = o'nuO'^'^, for the radial time-like geodesies 



'111/ 

comes out as 

2 mr) 
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where V{r) and Q{r) are are functions of r, given by 



V{r) = 2Kr^ {2C^G - l) + 2rWKr^ + 1 {G {G^G - l) + G^K) 
+r [AG^G - 3) + 2GWKr^ + 1 



Q(r) = (^Gr + Vi^r2 + 1 j ^ 



C2 I G'+^^5r!±l ) _i (34) 



The rotation parameter for the radial time-like geodesies, cu^ = w^j^cj^'^ = 0. 
Before we move on, let us make a few comments. First of all, note that O and 
diverge at r = 0, corresponding to the singularity of the BST at that point. These 
also diverges at (^\/K + + G) = 1, the turning point for outgoing radial 
time-like geodesies (see discussion after eq.( l25|) ) and indicates that the geodesies 
focus or de-focus at the turning point. 

Next we present the ESR parameters for the circular time-like geodesies. For 
these, the components are given as 

(G + 2GKr^ + 2KrVTTK^) 

B r = 

B\ = 
B\ = 



2r(l + Kr^fl^il + 2Kr^ + 2Gr^l + Kr'^f/^ 



2r{Gr + Vl + Kr^)il + 2Kr^ + 2GrVl + Kr^/^ 



[1 + 2Kr^ + 2GrVl + Kr^y/'^ 



^ iG + 2GKr' + 2KrVm^) 

rVl + Kr^{l + 2Kr^ + 2Gr^\ + Kr'^fl'^ ' 
and the diagonal elements of B^^ are all zero. This implies the expansion coeffi- 
cient G = for the circular time-like geodesies. The shear coefficient is 



a = : : : — . (OO) 

4rVl + Kr^{Gr + Vl + Kr^){l + 2Kr^ + 2Gr^\ + Kr'^Y 
The rotation parameter uP' for the circular time-like geodesies is 

4rVl + Kr\\ + 2Kr'^ + 2Gr^\ + Kr'^y 
To discuss the behavior of these parameters, let us first consider the case = 0. 
In this case, the shear and rotation parameters reduce to 

n^Gv^xf 



a 



4r2(Gr + l)(2Gr + 1) 



2 



4r2(2Gr + l)2 ^ 
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These are positive everywhere, for any value of /3 and go to zero in the hmit of 
infinite r (remember that G is positive). For G > and K < the analysis of the 
shear and the rotation parameter cannot be simply guessed by the expressions 
in eqs. (|36l) and ( |371) as in this case one can have an upper bound of the radial 
distance up to which circular time-like geodesies can be obtained in Type II BST. 

Finally, let us briefly discuss null geodesies. For the radial null geodesies the 
evolution tensor B^i, is defined as 

B^, = P\B^.P\ , (39) 

where the projection tensor Pap is 

Pali = 9ap + n^up + npUa ■ (40) 

Here Ua is the first integral of null geodesic equations i.e., u^V pUa = and 
UaU" = 0. The vector Ua satisfies the following conditions: 

n^n" = , n„M" = -1 , M^V/jn^ = . (41) 

Now the evolution tensor (in an effectively one dimensional space) becomes 

Bf,u = QP^u + + Cj^^ , (42) 

where 

= ^{B^,-B,^). (43) 
The vector n" is tangent to the inward radial geodesies. If we take 

= -1- 
2C" 



r 

n 



^pcJ{Kr^ + 1) [Gr + \/Kr^ + l) 



2C2 [Gr + V^r2 + l) 
n"^ = 0, (44) 

then these components satisfy all the conditions in Eq. ( HTj) . Using these we 
can find all the components of P^^. The only non zero component of P^^ turns 
out to be P*^^ which is given as 

P\ = l- (45) 
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Using the above component of the projection tensor one finds that there is only 



The above expression of the expansion variable for the null radial geodesic con- 
gruence for Type II BST shows that G has a singularity at r = where BST 
of Type II itself is singular. Unlike the outgoing time-like radial geodesies the 
expansion variable for the outgoing null radial geodesies do not have any other 
singularity. 

3.3 Analysis of BST Type II spacetime in terms of the 
ESR parameters 

Having obtained the ESR parameters for BSTs of type II, let us analyze them in 
some detail. First, we focus on radial time-like geodesies. In this case, the ESR 
parameters can be obtained by setting = in eq.( l32|) and in eqs. (l33|) . As 
mentioned in the discussion after eq.(|6]), we can consider two cases here, namely 
/3 = 1 and /3 < 1. Let us focus on the former case. Here, the expression for the 
expansion parameters assumes a simple form 



with cr^ given from eqs.( l33l) and ( IMl) . with K = 0. 

We wish to understand the behavior of O and cr as a function of the affine 
parameter, A. (Note, however that in our method, no initial condition on the 
parameter or a is possible, unlike the case where one integrates the full set 
of Raychaudhuri equations with given initial conditions [9], [10] (see also |llj). 
We will proceed, keeping this in mind). To this end, we first numerically solve 
eq. fl25|) (with chosen upper and lower limits of the affine parameter) to express the 
radial coordinate r as a function of A. This is then fed back in eq. fH7|) to obtain 
the behavior of ^ as a function of A. For illustration purpose, we have chosen 
G = 10~^, and set C = 1, and the lower and upper limits of the affine parameter 
A have been set to —0.2 and 1.5. The upper limit of A is chosen so that r varies 
from zero to the turning point of u"^, which can be seen to be r ~ 1 in this case 



one non-zero element of B'^^, namely, -B\^. Consequently, for the radial null 
geodesies, we obtain 





(47) 
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Figure 1: Numerical solutions for r (solid 
blue), (dotted red), and cr = ^J'o^i^Jo^ 
(dashed magenta) for radial time-like 
geodesies in type II BST, function 
of \ ioY K = 0. For details, see text. 



Figure 2: Numerical solutions for r (solid 
blue) and 9 (dotted red) for radial null 
geodesies in type II BST, function of 
A for the value K = —10^®. For details, 
see text. 



(these numbers are simply for illustration). The result is shown in fig.([T]), where 
the solid blue, dotted red and the dashed magenta curves correspond to numerical 
solutions for r, B and a respectively, as a function of the affine parameter, with 
the chosen initial condition. From the figure, we see that ^ is always negative, 
confirming the focusing theorem [6], [7] given the validity of the SEC of eq. (fTT]) . 
Also, diverges at the upper and lower limits of r, signaling the turning points 
of the geodesies. We find that the case 7^ 7^ for radial time-like geodesies follow 
the same qualitative behavior. 

To illustrate the case of null geodesies, we have taken K ^ 0. We have followed 
a numerical procedure similar to that alluded to above, and chosen K = — 10~^, 
G = 10~^, and C and /3 has been set to unity. Here, the lower and upper limits of 
the affine parameter has been set to —0.2 and 1.389 respectively. Using the same 
numerical procedure as above, we solve for the expansion parameter of eq. (H6|) . 
and this is illustrated in fig. ([2]), where we have multiplied by a factor of 10^ 
to display the curves on the same graph. Note that in this case the expansion 
parameter diverges at r = as expected from eq. (H6l) . 

Before we conclude, we briefiy mention the expansion parameter for the BST 
of type II in the four dimensional case. The details are unimportant here, and 
these can be worked out exactly like the case 6* = |. For radial time- like geodesies, 
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the expansion parameter is G = PA{r)/B{r) where 



Air) 



(ArWKr^ + 1 {C^G^ + G^K - G) + AKr^ {2G^G - l) 



+ 



r {8G^G - 5) + AGWKr^ + 1 




VKr^ + 1 



r 



) 



- 1 



(48) 



For radial null geodesies in four dimensions, the expansion parameter turns out 
to be twice that of eg. (1461) . The shear parameter for the radial time-like case 
yields a lengthy expression analogous to eqs. (l33|) and which we omit for 

brevity, and is zero for the radial null before. 

4 Conclusions and Discussions 

In this paper, we have considered the energy conditions and geodesic deformations 
in Bertrand space-times. For simplicity, we have chosen the BST of type II 
(eq.(l3])), although we believe that the qualitative results will remain unchanged 
even for a type I BST. As the metric of BST of type II contains three parameters 
we had to check the probable ranges of these parameters which gives the spacetime 
physical significance. We have explicitly checked the weak and strong energy 
conditions for type II BSTs, and verified their validity within certain ranges 
of parameters. The ESR parameters for both time-like and null geodesies in 
the BST type II spacetime are calculated as functions of the radial coordinate. 
For simplicity most of the calculations are done in the equatorial plane. We 
have not explicitly solved the Raychaudhuri equation in the above mentioned 
spacetime but indirectly obtained its solution by expressing the radial coordinate 
in terms of the affine parameter in the ESR variables. In such a situation the 
ESR variables becomes functions of the affine parameter and these variables are 
now valid solutions of the Raychaudhuri equation for geodesic deformations in 
type II BSTs. While analyzing the solutions of the Raychaudhuri equation for 
geodesic deformations in type II BSTs we have confirmed the focusing theorem 
numerically. 

The behavior of the solution of the Raychaudhuri equation for geodesic de- 
formations in type II BST for the time-like geodesic congruence is analyzed in 
this article with some suitable choice of parameter values K, G and /3. In the 
specific case chosen it is seen that generally the outgoing radial time-like/null 
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geodesies diverge from r — 0, which is a singular point in this spacetime where 
the Ricci scalar diverges. The outgoing radial time-like geodesies do not diverge 
indefinitely, they do converge again at some other value of r which turns out to 
be a turning point in type II BST. Where as an outgoing radial null geodesic 
congruence do not converge at any finite r but it does show focusing property. 
At the turning point for the time-like radial geodesies, spacetime is not singular 
but the radial component of the tangent vector to the radial geodesies vanish at 
this point. For the radial geodesies the rotation parameter is always zero but 
the time-like radial geodesies do show extreme shear near r = and the turning 
point. The circular time-like geodesies do not show any focusing behavior, as 
expected. But the circular orbits do show rotation. 

Our analysis points to the fact that BSTs can be thought of as interesting 
realistic examples of static, spherically symmetric space-times, which are asymp- 
totically non-flat, and allow for stable, periodic orbits. This might be signiflcant 
in astrophysical scenarios : for example, one might ask if a realistic space-time 
near a compact object can be modeled via BST of type II or I. In this article the 
main attention was given to the geodesies of BST of type II and its ESR variables 
to understand the effect of spacetime curvature and probable singularities. To 
properly utilize BSTs, one must also have to think of the source of such kind 
of spacetime's and in future works one needs to look at this important aspect. 
For completeness, it will be of interest to study, in the same manner, BSTs of 
type I and analyze the parameter space of this theory. We leave this for a future 
publication. 

Acknowledgements 

It is a pleasure to thank Sayan Kar for several useful correspondences. 

References 

[1] V. Perlick, Class. Quantum Grav., 9 (1992) 1009. 
[2] J. Bertrand, Compt. Rend. 77 (1873) 849. 

[3] H. Goldstein, "Classical Mechanics", 2°*^ edition, Narosa Publishing House, 
(1993) 

[4] R Kumar and K. Bhattacharya, Eur. J. Rhys. 32, 895 (2011) ArXiv: 
1103.3338 [physics . class-ph] . 



14 



[5] A. Ballesteros, A. Enciso, F. J. Herranz and O. Ragnisco, Class. Quan- 
tum Grav., 25, (2008) 165005. 

[6] R. M. Wald, "General Relativity," University of Chicago (2006). 

[7] E. Poisson, "A Relativist's Toolkit : The Mathematics of Black Hole Me- 
chanics," Cambridge (2004). 

[8] The literature on the subject is vast. For an excellent recent review, see 



S. Kar and S. SenGupta, Pramana 69 (2007) 49 ArXiv: gr-qc/0611123 



[9] A. Dasgupta, H. Nandan and S. Kar, Phys. Rev. D 79 (2009) 124004 
iarXiv: 0809 .30741 [gr-qc] . 

[10] R. Koley, S. Pal and S. Kar, Am. J. Phys. 71 (2003) 1037, ArXiv: 
|gr-qc/0302065| 

[11] S. Ghosh, A. Dasgupta and S. Kar, Phys. Rev. D 83, 084001 (2011) 
larXiv: 1008 .50081 [gr-qc]. 



15 



